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Abstract 

We consider one-dimensional quantum spin chain, which is called XX model 
(XXO model or isotropic XY model) in a transverse magnetic field. We study 
the model on the infinite lattice at zero temperature. We are interested in the 
entropy of a subsystem [a block of L neighboring spins] it describes entanglement 
of the block with the rest of the ground state. For large blocks entropy scales 
logarithmically jlUl II 1 j . We prove the logarithmic formula for the leading term 
and calculate the next term. We discovered that the dependence on the magnetic 
field interacting with spins is very simple: the magnetic field effectively reduce the 
size of the subsystem. We also calculate entropy of a subsystem of a small size. 
We also evaluated Renyi and Tsallis entropies of the subsystem. We represented 
the entropy in terms of a Toeplitz determinant and calculated the asymptotic 
analytically. 
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1 Introduction 



In this paper we study entropy of a spin chain. The entropy is the main object in 
thermodynamics and statistical physics. It is also interesting for information theory 
[H 121 EH- We study the how entropy of a subsystem scales with the size of the subsystem. 
We discover an interesting dependence on the magnetic field interacting with spins. The 
physical system we consider is XX model in a transverse magnetic field. Hamiltonian 
for this model can be written as 

N 

H xx (h) = - + aX+i + ha z n ). (1) 

n=l 

Here °n are Pauli matrix, which describe spin operators on n-th lattice site, 

h is the magnetic field and iV is the number of total lattice sites of spin chain (we 
take N — > oo in this paper). This model has been solved by E. Lieb, T. Schultz and 
D. Mattis in zero-magnetic field case [3] and by E. Barouch and B.M. McCoy in the 
presence of a constant magnetic field [3]. Some exact calculation of time-dependent 
properties also exists, for examples, see Ref. |E] by E. Barouch, B.M. McCoy and 
Refs. [Zj by M. Dresden and D.B. Abraham, E. Barouch, G. Gallavotti and A. Martin- 
Lof. We shall consider the model in thermodynamic limit N oo. The ground state 
is ferromagnetic for \h\ > 2 while it's critical for \h\ < 2. The ground state \GS) 
is unique. So the entropy of the whole infinite ground state is zero, but it can be 
positive for a subsystem [a part of the ground state] . We shall calculate the entropy of 
L neighboring spins. We shall call the first L neighboring spins as sub-system A and 
the rest as sub-system B. We shall consider von Neumann entropy (S(pa)), Renyi [2] 
and Tsallis [5J entropies for sub-system A: 

S von Neumann = S (Pa) = ~Tr{p A In p A ), (2) 

^Renyi = S<*(pa) = ^_ lnTr (PA)> a^l and a > 0. (3) 

Trp Q - 1 

^Tsallis = x _ a ( 4 ) 
Von Neumann entropy is the standard one. Tsallis and Renyi entropy may also be 
important for both information theory and statistical physics, see jS] and When 
a — > 1, Tsallis and Renyi entropy turns into von Neumann entropy. 1 Here p A is the 
1 Tsallis entropy and Renyi entropy are algebraically related: ^Xsallis = ~ — i-a _1 
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reduced density matrix of sub-system A: 

p A = Tr B (p AB ) 
and the density matrix of the whole system is 

PAB = \GS)(GS\ 

for zero temperature. Since calculations for von Neumann entropy and Renyi entropy 
are much similar, we give the detail calculation for von Neumann entropy only. The 
explicit result for Renyi entropy will be given without derivation. 

G. Vidal, J.I. Latorre, E. Rico, and A. Kitaev emphasized the role of the entropy of 
the subsystem in information theory [it describes entanglement of the subsystem with 
the rest of the ground state] . They showed that for subsystems of the large size L von 
Neumann entropy scales logarithmically S(pa) ~ (1/3) In L, see [121 HI]. In this paper 
we prove this formula and calculate the next term of asymptotic decomposition. We 
also evaluated Renyi entropy of the subsystem. Before we give the full derivation in 
the following sections, we first summarize our results here. We discovered that one can 
introduce a universal scaling variable: 



£ = 2L, 



'hV 



We consider the magnetic field to be less then critical value for \h\ < h c . When 
the magnetic field is larger then the critical value h c = 2 then the ground state is 
ferromagnetic [all spins are parallel], it does not have any entropy. We calculated von 
Neumann entropy and Renyi entropy of block spins for the magnetic field smaller then 
critical value. Here we present the expression for entropies for large £ and small £ 

+ ( 1 -£D 



S a (pA) 




if < £ < 1 

(a = 1) (5) 



,. ln£ + Tp if£»l 

Here is a constant defined in Eq. EH When a = 1, Renyi entropy S a (pA) becomes 
von Neumann entropy, the coefficient for log£ in large £ expression becomes | and 
Yj ^ becomes 



-t 



Ti = " dt — + 



< ; i «osh(t/2) 1 



3t ' tsinh 2 (t/2) 2sinh 3 (t/2) J 
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Following Ref. 0], let us introduce two Majorana operators 

1-1 1-1 
c 2 «-i = (II<K and <*u = 010*1 



(6) 



n=l 



n=l 



on each site of the spin chain. Operators c n are hermitian and obey the anti-commutation 
relations {c m , c n } = 2S mn . In terms of operators c n , Hamiltonian Hxx can be rewritten 

as 



N 



(C2nC2n+l — C2„_lC2„+2 + hC2 n -\C2n) ■ 



(7) 



n=l 



Here different boundary effects can be ignored because we are only interested in cases 
with N oo. This Hamiltonian can be subsequently diagonalized by linearly trans- 
forming the operators c n . It has been obtained |UE] (also see [TUJ [TT] ) that 

(GS\c m \GS) = 0, (GS\c m c n \GS) = 8 mn + t(B N ) mn . (8) 

Here matrix B m can be written in a block form as 



B 



N 



n n_j ... 
n x n ; 



H 



and n t = ±- r de e - iW g(9), 
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(9) 



where both and Q{9) (for iV — > oo) are 2x2 matrix, 



9(9) 
-9(0) 



9(0) 



(10) 



1, -k F < 9 < k F , 
-1, k F < 9 < (2tt - kp) 
and kp = arccos(|/i|/2). Other correlations such as (GS\c m ■ ■ ■ c n \GS) are obtainable by 
Wick theorem. The Hilbert space of sub-system A can be spanned by 11^=1 {^j^ 1 0)^, 
where af is Pauli matrix, takes value or 1, and vector \Q) F denotes the ferromag- 
netic state with all spins up. We are also able to construct a set of fermionic operators 
bi and bf by defining 



2L 



J2 v ™nCm m = 1, • • • , 2L; b t = (d 2 i + ida+i) / 2 > ^ = 1 r--,L (11) 

n=l 

with v mn = (V) mn . Here matrix V is an orthogonal matrix. It's easy to verify that d m 
is hermitian operator and 

bf = (d 2l -id 2l+1 )/2, {b i ,b j } = 0, {&+&+} = 0, {!>,:. b,} <),. r (12) 
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In terms of fermionic operators bi and bf, the Hilbert space can also be spanned by 
n^=i{^i"} Pi |0)uac- Here pj takes value or 1, 2L fermionic operators 6,, bf and vacuum 
state |0}„ ac can be constructed by requiring 

k\0) vac = 0, / = 1,-..,L. (13) 

We shall choose a specific orthogonal matrix V later. 



2 Density Matrix of Sub-system A 

Let {ipi} be a set of orthogonal basis for Hilbert space of any physical system. Then 
the most general form for density matrix of this physical system can be written as 

p = J2c(i,jMi)(^j\. (i4) 

I, J 

Here c(J, J) are complex coefficients. We can introduce a set of operators P(I, J) by 
P(I, J) oc \ipi)(ipj\ and P(I, J) satisfying 

P(I,J)P(J,K) = 6 ItK \1> J )(il> I \, P(I,J)P{J,K) = 6 I)K \>t/> I )(1> J \. (15) 

There is no summation over repeated index in these formula. We shall use an explicit 
summation symbol through the whole paper. Then we can write the density matrix as 

p = Y,c(I,J)P(I,J), c(I,J)=Tr(pP(J,I)). (16) 
i,j 

Now let us consider quantum spin chain defined in Eq. ^ For the sub-system A, the 
complete set of operators P(I, J) can be generated by Yli=i 0%- Here we take operator 
Oi to be any one of the four operators {bf ,bi,bfbi,bibf} (Remember that bi and bf 
are fermionic operators defined in Eq. ITTj) . It's easy to find that P(J, I) = (Il^=i Orf 
if P(I, J) = Il£=i Oi. Here f means hermitian conjugation. Therefore, the reduced 
density matrix for sub-system A can be represented as 

p A = Y,Tr AB (pab(1[oA no,. (17) 
V i=i / i=i 

Here the summation is over all possible different terms E[i=i Oi. For the whole system 
to be in pure state \GS), the density matrix pab is represented by \GS)(GS\. Then 
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we have the expression for p A as following 

p A = Y,(GS\(f[O i y\GS)l[O l . (18) 

i=l i=l 

This is the expression of density matrix with the coefficients related to multi-point cor- 
relation functions. These correlation functions are well studied in the physics literature 
Now let us choose matrix V in Eq. El so that the set of fermionic basis {bf} and 
{bi} satisfy an equation 

(GSlbibjlGS) = 0, (GSlbfbjlGS) = 5 hj (GS\btbi\GS). (19) 

Then the reduced density matrix p A represented as sum of products in Eq. [TH] can be 
represented as a product of sums 

L 

Pa = l[((GS\bth t \GS)bth t + (GSlhbtlGS^bt). (20) 
i=i 

Here we used the equations (GS\bi\GS) = = (GS\bf\GS) and Wick theorem. This 
fermionic basis was suggested by G. Vidal, J.I. Latorre, E. Rico and A. Kitaev in 
Ref . fOJ ■ A similar result for the density matrix of a subsystem in terms of free 
spinless fermion model was obtained by C.A. Cheong and CL. Henley in Ref. [T2] . 



3 Closed Form for The Entropy 



Now let us find a matrix V in Eq. ^2 which will block-diagonalize correlation functions 
of Majorana operators c n . From Eqs. ^2 and EH we have the following expression for 
correlation function of d n operators: 



2L 2L 



(GS\d m d n \GS) = J^J^VmiiGSlciCjlGS}-* 

i=l j=l 

{GS\c m c n \GS) = 5 mn + i(Bij mn , 
(GS\d m d n \GS) = 5 mn + i(Bi J ) ni „. 



(21) 



The last equation is the definition of a matrix Bl. Matrix Bl is the sub-matrix of Bn 
defined in Eq. El with m, n — 1, 2, . . . , L. We also require Bl to be the form ^UJ 



B L = VB L V q 





f 






< 










= 




:) 




V" 1 






I- 1 





(22) 



Here matrix O is a diagonal matrix with elements v m (all v m are real numbers). There- 
fore, choosing matrix V satisfying Eq. |22in Eq. ^2 we obtain 2L operators \b{\ and 
{bf} with following expectation values 

1 + u m 



(GS\b m \GS) = , (GS\b m b n \GS) = , (GS\b+b n \GS) = 5„ 
Using the simple expression for reduced density matrix pA in Eq. EDI we obtain 

pa = n C-^btk + l —^hbi 



(23) 



1=1 



(24) 



This form immediately gives us all the eigenvalues X Xl x 2 ---x L of reduced density matrix 
PA, 



A 



x 1 x 2 --xi J 



n 

i=l 



i + -ir^ 



, Xj = 0, 1 Vz. 



(25) 



Note that in total we have 2 L eigenvalues. Hence, the entropy of pa from Eq. |21becomes 

S(pa) = I>(l,z/ m ) (26) 

with 



m=l 



. . X+V X + v x—v x—v 
= — ln(— ^— ) — ln(^— ). 



(27) 



We shall use this result further to obtain analytical asymptotic. Function e(l, v) in 
Eq. |2E1 is equal to the Shannon entropy function H(^j^). However, in the following 
calculation (Eq. I32|). we will need the more general function e(x,u) instead of H{y). 
Notice further that matrix Bl can have a direct product form, i.e. 



B L = G L (8) 



1 
-1 



with Gl 



go g-i 
gi go 



9l-L 



V 9L-1 



(28) 



where g\ is defined as 



9l = ±- r d6e- iW g(6) and g{6) 

Z7T JO 



1, —k F <9 < hp, 
•1, k F < 9 < (2tt - k F ) 



(29) 
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and kp = arccos(|/i|/2). From Eqs. l2*2l and 1281 we conclude that all just the 

eigenvalues of real symmetric matrix Gl. 

However, to obtain all eigenvalues v m directly from matrix Gl is a non-trivial task. 
Let us introduce 

D L (A) = det(G L (A) = AJ L - Gl) . (30) 

Here Gl is a Toeplitz matrix (see JTj) and Jl is the identity matrix of dimension L. 
Obviously we also have 

L 

£>l(A) = l[(\-Vm). (31) 
m=l 

From the Cauchy residue theorem and analytical property of e(x,v), then S(pa) can 
be rewritten as 

S(p A ) = Hm lim — I e(l + e, A)d In D L (A) . (32) 

e-^0+ 5^0+ ZTC1 Jc(e,5) 

Here the contour c(e, 5) in Fig [T] encircles all zeros of D^(X), but function e(l + e, A) is 
analytic within the contour. Just like Toeplitz matrix Gl is generated by function g{9) 
in Eqs. EH] and EH] [see next section], Toeplitz matrix Gl(A) is generated by function 
g{9) defined by 

!A - 1, -k F < 9 < k F , 
(33) 
A + 1, k F < 9 < (2tt - k F ). 

Notice that g(8) is a piecewise constant function of 9 on the unit circle, with jumps at 
9 = ±k F . Hence, if one can obtain the determinant of this Toeplitz matrix analytically, 
one will be able to get a closed analytical result for S(pa) which is our new result. 
Now, the calculation of S(pa) reduces to the calculation of the determinant of Toeplitz 
matrix Gl(A). Before we calculate the determinant of Toeplitz matrix Gl(A), we also 
want to point out two special cases which allow us to obtain an explicit form for these 
eigenvalues v m and hence the entropy S(pa)- These are cases with small lattice size 
of subsystem A and magnetic h close to critical values ±2. Let us take k F = k F for 
k F < | or k F = vr — kp for k F > |. For cases kp L 1, Toeplitz matrix Gl can be well 
approximated by a matrix with diagonal elements (2kp/iT — 1) and all other matrix 
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Figure 1: The contour c(e, S). Bold lines (—00, — 1 — e) and (1+e, 00) are the cuts of integrand 
e(l + e, A). Zeros of -Dl(-X) (Ea.l5T|) are located on bold line (—1, 1) and this line becomes the 
cut of dlog -Dl(A) for L — > 00 (Eq. I48j) . The arrow is the direction of the route of integral we 
take and R is the radius of circles. 

elements equal to 2kp/ir. Hence, we can obtain all eigenvalues for Toeplitz matrix Gl 
as {2Lkp/TT — 1, —1, —1, • • • , —1} and S(pa) becomes 

S(p A ) « In -^U, < k F h < 1. (34) 

7T 2L/c F 

Expression above can also be re-expressed in terms of h as 

S(p A )« 2L(1 - fe2/4) * ln \ < (1 - /> 2 /#L « 1. (35) 

7T 2L(1 - /i 2 /4)2 

4 Toeplitz Matrix and Fisher-Hartwig Conjecture 

Toeplitz matrix T L [0] is said to be generated by function (f)(9) if 

T L [0] = i,j = l,---,L-l (36) 

where 

j ^e^dO (37) 



Z7T JO 
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is the Z-th Fourier coefficient of generating function (f)(9). The determinant of 2l[</>] 
is denoted by D^. The asymptotic behavior of the determinant of Toeplitz matrix 
with singular generating function was initiated by T.T. Wu ^13] in his study of spin 
correlation in two-dimensional Ising model and later incorporated into a more general 
conjecture, i.e., the famous Fisher-Hartwig conjecture [T4l Ho^ ITol ITTj. 
Fisher-Hartwig Conjecture: Suppose the generating function of Toeplitz matrix 
(f)(8) is singular in the following form 

m = mf[tp r ,8Mu ar , 9r (9) (38) 

r=l 

where 

tpr,eM = exp[-i/3 r (7T-0 + r )], 6 r < 9 < 2vr + 6 r (39) 
u ar ,eM = (2-2cos(0-0 r )) ar , ^a T >- 1 - (40) 

and i/j: T — > C is a smooth non-vanishing function with zero winding number. Then 
as n — > oo, the determinant of T L [0] 



^ L = WD L fn La?_/3? ) *[V>, W, {A}, {0,}], l ^ oo. 



(41) 



Here ^ r [V'] — exp (^/ o 27r ln^(0)d0). Further assuming that there exists Weiner-Hopf 
factorization 

if)(6) =^]V+(exp(i0))V-(exp(-i0)), (42) 
then constant £[0, {A}, {0«}] i n Eq. HTJcan be written as 

R 

£[iP, {«,},{A},{^}] = ^[^n G ( 1 +^ + A)G'(l + « J -A)/G'(l + 2« J ) 

i=l 

x n(V'-(exp(i0 i ))J ^+(exp(-i0O)J 

x [] l-exp(i(^-^)) , (43) 

G is the Barnes G-function, £[0] = exp(X)fcli kskS-k), and is the fc-th Fourier coef- 
ficient of In -0(0). The Barnes G-function is defined as 

oo 

G(l + z) = (2rr) z / 2 e^ z+l)z/2 -^ z2/2 J] {(1 + 2/n)" e - 2+22/(2n) }, (44) 

n=l 
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where 7^ is Euler constant and its numerical value is 0.5772156649 • • •. This conjecture 
has not been proven for general case. However, there are various special cases for which 
the conjecture was proven. 

For our case, the generating function g{9) has two jumps at 9 = ±kp and it has 
the following canonical factorization 

m = W)fam,kMt*w,-kM ( 45 ) 

with 

^) = (A + l)(^f I J , 0(A) = -ft(A) = ft(A) = — In^fy. (46) 

The function t was defined in Eq. EH1 We fix the branch of the logarithm in the following 

way 

- 7T < arg [j^jj < n ■ ^ 

For A i [-1,1], we know that |3?(/3i(A))| < 1 and |3fc(0 2 (A))| < \ and Fisher-Hartwig 
conjecture was PROVEN by E.L. Basor for this case Therefore, we will call it 
the theorem instead of conjecture for our application. Hence following the theorem in 
Eq. UU the determinant D L (\) of AIl — Gl can be asymptotically represented as 

D L (X) = (2-2cos(2A; F ))" /32(A) {G'(l+/3(A))G'(l-/3(A))} 2 

(A + 1)((A + 1)/(A - l)Y kF/n ] L L- 2 ^ x \ (48) 

Here L is the length of sub-system A and G is the Barnes G-function and 

G(l +0(A))G(1 -0(A)) = e-^+^W ft { (l - ■ 



n=l 



(49) 



5 Asymptotic Form of The Entropy 

Now, let us come back to the calculation of entropy S(pa). For later convenience, let 
us define 

r(A) -i;^8W' (50) 
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Taking logarithmic derivative of -Dl(A) (Eq. 148)1. we obtain 
dlnD L (A) fl-k F /n k F /n^ 



dX \ 1 + A 1-A 

4 P(X) 



IL 



j_ ( , + A)( , A) V lnL + ln ( 2 l sin ^D + C 1 + 7is) + T(A)). (51) 

Let us substitute the asymptotic form above for dlnDL(A)/dA into Eq. for entropy 
S(pa): 

S( PA ) = hm lim -1 / e(l + e, A) f 1 = - ^ L+ 
^ e ^o+ <5^o+ 2vri J c (e,5) V \ 1 + A 1 - Ay 

lim lim-^/ dA ^ 1 + ''^^ (lnL + ln(2| sin k F \) + (1 + 7js ) + T(A)) ; (52) 
e ^ + 5^0+ 7T 2 Jc(e,8) (1 + A)(1-A) V V ' 17 V ' v /y>v 7 

where the contour is taken as shown in Fig.^ The first integral in Eq.|32]can be carried 
out by using the residue theorem and the definition of function e(x, v) in Eq. EH We 
found that the linear term in L for entropy S(pa) vanishes. The second integral can 
be calculated as follows: First, we notice that 



c( «,, dA( "' ) = l/Af + /FEl5 + ii3a + /cB3j dA( '" ) (53) 

Second, we can show that the contribution of integral from the circular arcs FEl3 and 

CbA vanishes. Therefore, the entropy (Eq. [52} can be written as 

c, f n ,• 2 / /-i+io+ a+io- e (l + e,A)/?(A) 
S(p A ) = hm — / + / dA 



e^0+ 7T 2 Vl+iO+ i-n-io-' (1 + A)(1-A) 
x (lnL + ln(2|sinfc F |) + (l + 7£) + T(A)). (54) 



For further simplification, we shall use the fact that 



(3(x + i0 ± ) = — In T i(vr - 0+) = -iW(x) + (- - 0+) (55) 

2m \ 1 — x J 2 



for x G (—1, 1) and 



W(x) = ^\n^-?-. (56) 



We can now write the entropy S(pa) as 

S(p A ) = - 2 f dx^%4( lnL + ln (2|sinfc F |) + (l+ 7 i ; ; 

" 2TI- 1 /* e(l, x) ( (| + i^(x)) 3 (| - W{x) f \ 

J-i 1-x 2 \n 2 -{l + iW{x)f n 2 -(\-\W(x)) 2 ) k ' 
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where e(l,x) is defined in Eq. |23 This expression for S(pa) contains two integrals. 
The first integral can be done exactly as 

2/-1.1+X.1 + X 1-x.l-x, 1 1 , ro , 

— / ax { In In = -. 58 

W-i v 2 2 2 2 J l-x 2 3 V ; 

The second integral in Eq. |37| can be written as 

^ 7T 2 7-1 v l-x 2 1 + x 2 ; 
/ {\ + iW{x)f {\-iW{x)f \ 

[n 2 - (i + W(a;)) 2 n 2 - (| - iW(x)) 2 J ' 1 ' 

which can be further simplified ^H]- Finally we have that 



S{p A ) = \\nL + UnU-(^ j+^ + Ti, L - oo (60) 

with 

f» (e" 1 1 cosh(t/2) 1 

1 " Jo \"3T + tsinh 2 (i/2) ~ 2sinh 3 (t/2) J ' ( } 

for XX model. The leading term of asymptotic of the entropy |lnL in Eq. IBTH was 
first obtained based on numerical calculation and a simple conformal argument in 
Ref. fOm] m the context of entanglement. We also want to mention that a complete 
conformal derivation for this entropy was found in Ref. JH]- One can numerically 
evaluate Ti to very high accuracy to be 0.4950179 • • •. For zero magnetic field (h = 0) 
case, the constant term Ti + In 2/3 for S(pa) is close to but different from (7r/3)ln2, 
which can be found by taking numerical accuracy to be more than five digits. 



6 Summary 

In this paper, we study the entropy of a block of L neighboring spins in XX model 
with the presence of the transverse magnetic field. We obtain Eq. ESI and Eq. EH 
for the von Neumann entropy of a block of L neighboring spins in XX with small L 
and large L respectively. It's interesting to note that there is a natural length scale 
L;,, = 1/ — (j^J ^ for \h\ < 2 to incorporating the magnetic field effects. When \h\ 
increases from less than 2 into larger than 2, the system evolves from critical phase 
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into ferromagnetic phase. The ferromagnetic phase does not have any entropy. The 
scale Lh shows singular behavior at the critical value of the magnetic field h c = 2. We 
discovered that one can introduce the universal scaling variable £ = 2L/L/j: 



C = 2L 1 



for \h\ < 2. Then we can express the von Neumann entropy of L neighboring spins in 
following simple form: 

, £ In £ if < C < 1 

%)= : (62) 




llnC + T, if£»l 



with 



T = n&tl— 1 cosh(t/2) \ 

1 ' Jo \3t tsinh 2 (t/2) 2sinh 3 (t/2) J ' 

For small lattice and magnetic field close to ±2, we obtain the result directly. To 
obtain the result for large L asymptotically, we first expressed the entropy in terms of 
the determinant of a Toeplitz matrix. Then we used a special case of Fisher-Hartwig 
conjecture ^3] and this special case was PROVEN in [To] . 

From similar calculation, we also obtain the Renyi entropy in Eq. |3] to be 

, ^-lnfff ) a + (1 - f Y) ifO<£<l 

S*( P a) = VW s \ 2nJ > (63) 



i±^ln/: + TS Q} if£»l 



Here 



T<«> = -l| 1 i dxi^(«(i-i^)) + «(i + i^W)), (64) 

- r^'«((T)" + (V)")- < 65 > 

* W . iL hirw ._ 7E + 1 (_!___!_), (66) 



with 7£ the Euler Constant and T(x) the well-known Gamma Function. 



14 



Acknowledgments 

We would like to thank Prof. B.M. McCoy, Prof. A. Abanov and Mr. F. Franchini for 
useful discussions. This work was supported by National Science Foundation (USA) 
under Grant No. DMR-0073058 and PHY-9988566. 

References 

[1] C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schumacher, Phys. Rev. A 53, 
2046, (1996) 



[2 
[3 
[4 
[5 
[6] 
[7] 



[9 

[10 

[11 
[12 

[13 

[14 



A. Renyi, Probability Theory, North-Holland, Amsterdam, 1970 

S. Abe and A. K. Rajagopal, Phys. Rev. A 60, 3461, (1999) 

E. Lieb, T. Schultz and D. Mattis, Ann. Phys. 16, 407, (1961) 

E. Barouch and B.M. McCoy, Phys. Rev. A 3, 786, (1971) 

E. Barouch, B.M. McCoy and M. Dresden, Phys. Rev. A 2, 1075, (1970) 

D.B. Abraham, E. Barouch, G. Gallavotti and A. Martin-L6f, Phys. Rev. Lett. 
25, 1449, (1970); Studies in Appl. Math. 50, 121, (1971); ibid 51, 211, (1972) 

C. Tsallis, J. Stat. Phys. 52, 479, (1988) 

N.M. Bogoliubov, A.G. Izergin, and V.E. Korepin, Quantum Inverse Scattering 
Method and Correlation Functions, Cambridge Univ. Press, Cambridge, 1993 

G. Vidal, J.I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90, 227902, (2003) 

J.I. Latorre, E. Rico, and G. Vidal, arXiv: quant- ph/0304098 



C.A. Cheong and C.L. Henley, arXiv: jcond-mat/0206196 
T.T. Wu, Phys. Rev. 149, 380, (1966) 

M.E. Fisher and R.E. Hartwig, Adv. Chem. Phys. 15, 333, (1968) 



15 



[15] E.L. Basor, Indiana Math. J. 28, 975, (1979) 

[16] E.L. Basor and C.A. Tracy, Physica A 177, 167, (1991) 

[17] A. Bottcher and B. Silbermann, Analysis of Toeplitz Operators, Springer- Verlag, 
Berlin, 1990 

[18] F. Franchini, Private communication 

[19] V.E. Korepin, arXiv: |cond-mat/0311056| 



Appendix: Simplification of Formula 



In this appendix, we show more details for simplification of To (JHH|) in detail. In order 
to simplify T , we will use the Function ip(x), which is defined as 

d 



ip(x) = -j- In T(x) 
ax 



OO -y OO 1 



n + x 



(6? 



n=0 1 n=0 

with the Euler Constant and T(x) the well-known Gamma Function, and the prop- 
erty 



1 



ip(x + 1) = ip(x) H — . 



x 



Introducing z (z) 



[—) iW(x) and using Eqs. (jUH|) and (JB9J) . we obtain 



\ + iW(x] 



(± - iW(x) 



n=l 



- - 1 1 + iW(x)) 2 + n 2 - (| - iW(x) N 2 



2 

= ^(l)-l-^(i-ilV(x))-^(i + iiy(x) 
by using Eq. (jnHj) and definition for z and z. Hence, we obtain 



1 r 1 
T = — dx 



7T 2 J-l 



1 1 + X 1 1 — X 

in In 



Tx 



l-x 2 l+x 2 
^(1) - 1 - ~ - iW(x)) - + iW(x 

1 + lE 



T V2 



2 r V2 



(69) 



(70) 



(71) 
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with Ti defined as 
Ti = 



l r 1 



2tt 2 



dx 



1 1 + x 1 1 — x 
m In 



l-x 



x 



2 1 + x 
^ iW{x))+^- + iW{x)) 



^2 v ^2 
We now perform a change of variable using w = ^ In 

-2 



—2 r 00 

Ti = — / dw (In [2 cosh(7rw)] — nw tanh(7rw)) 

TV JO 

2 / r V2 



1 1 



We note that 



In [2 cosh(irw)] - ttw tanh(7r«;) = ^1 - In (l + e " 2 ™ Q ) 



a=l 



Hence we can rewrite 
-2i f°° 



-2i f°° ( d \ T ( \ - \w 

Ti = / dw (In [2 cosh(7ru;)l — mu tanh(7rw)) • - — In — 7 . 

Wo 11 1 ,] K " \dw) T{\ + iw) 

Using the following expression for the Logarithm of the Gamma Function: 



lnT(z 



Jo 



z-1 



1 - e"^" 1 )* 



^dt 
t 



1 - e-* 

which is particularly convenient because we need only the imaginary part of it: 



In 



r (I - ™) 

r(i + i 



Jo 



2we~ 



sin(iui) 
sinh(t/2) 



dt 



After some elementary but tedious calculation, finally we obtain 



Ti 



Jo 



dt 



+ 



COS 



Htm \ 



3t tsinh 2 (t/2) 2sinh 3 (t/2)j 
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